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Abstract 

After briefly reviewing (3/3 decay as a test of the neutrino mass, I examine 
s * | the nuclear structure involved in this process. Simple formulas (a la Pade) are 
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designed for the transition amplitudes and the general behavior of (3(3 decay 
amplitudes in the quasiparticle random phase approximation are discussed. 
Results of a calculation for 48 Ca, 76 Ge, 82 Se, 100 Mo, 128 Te and 130 Te nuclei are 
presented, in which the particle-particle interaction strengths have been fixed 
by invoking the partial restoration of the isospin and Wigner SU(4) symmetries. 
An upper limit of < m v >~ 1 eV is obtained for the effective neutrino mass. 

o 

g : PACS numbers: 23.40.Hc, 21.10.Re, 21.60.Jz 

1. Introduction 

The double beta {(3(3) decay is a nice example of the interrelation between the 
Particle Physics and the Nuclear Physics: we can get information on the properties of 
the neutrino and the weak interaction from the (3(3 decay only if we know who to deal 
we the nuclear structure involved in the process. There are already several well-known 
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reviews on the neutrino physics ]T|, H and therefore, after a brief historical overview, 
I will limit this talk to the nuclear facet of the problem. More precisely, some recent 
developments performed by our group will be summarized. 

Because of the pairing force there are approximately 50 nuclear systems in which 
the odd-odd isobar, within the isobaric triplet (A,Z), (A, Z + 1), (A,Z + 2), has 
a higher mass than its neighbors. Within such a scenario the single (3 decay, is 
energetically forbidden and the initial nucleus can disintegrate only via the (3(3 decay. 
This is a second-order weak interaction process, similar to electromagnetic processes 
such as the atomic Raman scattering and the nuclear 77 decay ||. 
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The modes by which f}(3 decay can take place are connected with the neutrino [i>) 
- antineutrino (u) distinction. If v and v are defined by the transitions: 

n — > p + e™ + v (1) 
v + n — > p + e~, 

the decay (A, Z) — > (A, Z + 2) can occur by successive (3 decays: 

(A,Z) — ► (A,Z + l) + e- + v 

— ► (A,Z + 2) + 2e~ + 2v, (2) 

passing through the intermediate virtual states of the (A, Z + 1) nucleus. 

Yet the neutrino is the only fermion in lacking a additionally conserved quantum 
number that differences between v and v. Thus it is possible that the neutrino is a 
Majorana particle, i.e., equal to its own antiparticle (a la 7T°). [] When v = v the 
process 

(A, Z) — > (A,Z + l) + e~ + u = {A,Z + 1) + e" + u 

— > (A,Z + 2) + 2e~ (3) 

is also allowed. In absence of the helicity suppression (as would be natural before the 
observation of parity violation) this neutrinoless (/3/3q v ) mode is favoured by phase 
space over the two-neutrino {fifi-iv) mode by a factor of 10 7 — 10 9 : T 0u ~ (10 13 — 10 15 ) 
y while T 2u ~ (10 20 — 10 24 ) y. Several searches for the /3/3 decay has been made by 
the early 1950s with the result that T > 10 17 . This seemed to point that v ^ v and 
prompted the introduction of the lepton number L to distinguish v from v: L — +1 
was attributed to e~ and v and L = — 1 to e + and v. The assumption that the 
additive lepton number is conserved then allows the /?/3 2l , decay but prohibits the 
/3/3o^ one, for which AL = 2. 

But with the discovery in 1957 that the parity is not conserved for the weak 
interaction it was realized that the Majorana/Dirac character of the neutrino was 
1 A Dirac particle can be viewed as a combination of two Majorana particles with equal mass and 
opposite CP properties and their contribution to the f3(3o v decay cancel. 
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still in question. If 

n — > p + e~ + Pre (4) 
vlh + n — > p + e~ 

then the second process in fl3|) is forbidden because the right handed neutrino has the 
wrong helicity to be reabsorbed. Therefore the 75- invariance of the weak interaction 
could account for no (3(3 0v decay, regardless of the Dirac or Majorana nature of the 
neutrino. Otherwise, this decay can be observed only when the lepton number is not 
conserved and the neutrino is a massive Majorana particle. This event discouraged 
experimental searches for a long time, but with the development of modern gauge 
theories the situation began to change. In fact, there are many reasons for the re- 
naissance of interest in /3/9-decay over the past decade. The most important one is 
that, if there is any new physics beyond the standard SU(2)l x U(l) gauge model of 
electroweak interactions, the j3j3 0u decay will play a crucial role in shaping the ulti- 
mate theory. Moreover, no solid theoretical principle precludes neutrinos from having 
mass and the most attractive extensions of the standard model require neutrinos to 
be massive. The theory is neither capable of predicting the scale of neutrino masses 
any better than it can fix the masses of quarks and charged leptons. 
Since the half-lives can be cast in the simple form 

(where Q's and Ai's are, respectively, the phase space factors and nuclear matrix 
elements and < m u > is the effective neutrino mass) it is clear then that we shall not 
understand the (3{3q u decay unless we understand the Pfov decay. The last is one of the 
slowest process observed so far in nature and offers a unique opportunity for testing 
the nuclear physics techniques for half-lives > 10 20 y. Thus, the comprehension of the 
(3(3 transition mechanism cannot but help advance knowledge of physics in general. 

It is worth noting that more than 30 (3(3 decay experiments are underway or in 
stages of planning and construction. Until now positive evidence of the /3/3 2i , decay 
2 We assume for simplicity that weak interactions with right-handed currents do not play an 
essential role in the neutrinoless mode. 
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mode has been found for 76 Ge, 82 Se, 100 Mo, 128 Te, 130 Te and 238 U. Yet, despite the 
colossal experimental progress the neutrinoless, lepton violating decay, if it exists, has 
escaped detection until now. [] 

The (3(3 decays occur in medium-mass nuclei that are rather far from closed shells, 
and we all know that shell-model calculations are practical only when the number 
of valence nucleons is relatively small. Therefore, at the present time, the nuclear 
structure method most widely used is the quasiparticle random phase approximation 
(QRPA). Within this model the /5/?-decay amplitudes are very sensitive to the inter- 
action parameter in the particle-particle (PP) channel, usually denoted by g pp . It is 
still more interesting that, close to the expected value for g pp , the (3(3 matrix elements 
go to zero. But when a physical quantity has a zero (or near zero) a conservative law 
should, very likely, be at its origin. Thus, resorting to a toy model, I will first discuss 
the general behavior of the nuclear matrix elements || [(J. Later, I will show that 
they have zeros because of the restoration of the isospin and Wigner S77(4) symmetry 
H [Tj. This is not surprising since the Fermi (F) and Gamow- Teller (GT) operators r± 
and cn±, relevant in the (3(3 decay, are infinitesimal generators of SU(2) and S77(4), 
respectively. Finally, we use the concept of restoration of these symmetries to fix the 
PP interaction strengths and to estimate the (3(3 matrix elements. It can be argued 
that, the 577(4) symmetry is badly broken in heavier nuclei and that therefore our 
recipe is quite arbitrary. I will show, yet, that the residual interaction is capable to 
overcome most of the 577(4) breaking caused by the spin-orbit splitting. 

2. General behavior of M. 2 v and _M „ in the QRPA 

Independently of the nucleus that decays, of the residual interaction that is used, 
and of the configuration space that is employed, the /3/3-moments as a function of g pp 
always exhibit the following features: 

(i) The 2v moments have first a zero and latter on a pole at which the QRPA 
collapses. 

(ii) The zeros and poles of M.q v for the virtual states with spin and parity J 71 = 1 + 
are strongly correlated with the zeros and poles of M.i v . 

3 The P/3 0l , sensitivities have changed from ~5x 10 15 in 1948 to ~ 5 x 10 23 in 1987. 
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(iii) The total (5(3q v moments also possess zeros but at significantly larger values 
of g pp . 

The behaviour of the (3(3 moments for several nuclei are illustrated in Fig. [T| 
These results have been obtained with a 5 force, using standard parametrization [[/J. 
Instead of the parameter g pp , I use here the ratio between the triplet and singlet 
coupling strengths in the PP channel, i.e., t = vf p /v pp . Calculations with finite range 
interactions yield similar results ||. 

I will resort now to the single mode model (SMM) description []|, || of the (3(3- 
decays in the 48 Ca — > A8 Ti and 100 Mo — > 100 Ru systems. This is the simplest version 
of the QRPA with only one intermediate state for each J n . It allows to express the 
moments a la Alaga i.e., as the unperturbed matrix elements hA\ v and Ai^ u (J + ) 
multiplied by the effective charges: 

M „{J+) = Ml(J + )— (l+^p), (6) 

UJJ+ \ UJ" J 

Here G(J + ) = G(pn,pn; J + ) are the PP matrix elements (proportional to t (or to 
g pp )), u>° is the unperturbed energy, and ojj+ are the perturbed energies. It will be 
assumed here that the isospin symmetry is strictly conserved, in which case (as it will 
be seen latter on) A / i2^(0 + ) = M.q u {Q + ) = 0. When the pairing factors are estimated 
in the usual manner, one gets 

io = oj°Jl + F(34 + 9F/u°)/25(jj° + 16G(1 + F /uj°)/25uj , (7) 

for the single pair configurations [0f7/2(n)0f 7 / 2 (p)]j+ in 48 Ca and 

u = uj°^/i + 4F(45 + F/uj°)/225uj + G(270 + 172F/lu° + 49G/cj )/225w , (8) 

for [0^7/2 (ti) 0^9/2 (p)]j+ in 100 Mo. Therefore, while the numerators in Eqs. @ and @ 
depend only on the PP matrix elements, their denominators depend on the particle- 
hole (PH) matrix elements F(J + ) = F(pn,pn; J + ), as well. The numbers in the last 
two equations arise from the pairing factors. As illustrated in Fig. 0, the SMM is a 
fair first-order approximation for the ftfov decays in 48 Ca and 100 Mo nuclei. 
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The role played by the ground state correlations (GSC) in building up Eqs. @ 
and @ can be summarized as follows: 

(a) The numerator, i.e., the factor (1 + G/u°), comes from the interference between 
the forward and backward going contributions. These contribute coherently in the 
PP channel and totally out of phase in the PH channel. 

(b) The G 2 and F 2 terms in the denominator are very strongly quenched by the GSC, 
while the GF term is enhanced by the same effect. In particular, for 48 Ca the term 
quadratic in G does not contribute at all. 

It can be stated therefore that, within the SMM and because of the GSC, the M.2u 
matrix element is mainly a bilinear function of G(l + ). Besides, it passes through 
zero at G(l + ) = — uj° and has a pole when uj\+ = 0. Similarly, all AA.q v {J + ) moments 
turn out to be quotients of a linear function of G(J + ) and the square root of another 
linear function of G(J + ). Both the zero and the pole of }Aq v {1 + ) matrix element 
coincide with those of the 2v moment. Besides, as the magnitudes of G(J) and F(J) 
decrease fairly rapidly with J (see Table [l]), the quenching effect, induced by the PP 
interaction, mainly concerns the allowed Ov moment. For higher order multipoles it 
could be reasonable to expand the denominator in Eq. in powers of G(J + )/u 
and to keep only the linear term. This term strongly cancels with a similar term in 
the numerator and the net result is a weak linear dependence of the M.$ V (J + 7^ 1 + ) 
moments on the PP strength. Obviously, for the last approximation to be valid, the 
parameter t (or g pp ) has to be small enough to keep 0;!+ real. Briefly, the SMM can 
account for all four points raised above, and leads to the following approximations 

M 2v ^M 2 „(t = 0)\^^, (9) 
1 — tjt\ 



and 



M u = M 0u (J* = l + ;t = 0) ' 1/1,1 



V 1 " t/ti 

+ M 0v (r ^l + ;t = 0)(l-t/t 2 ), (10) 

where t\ > t and ti 3> tj., and the condition t < t\ is fulfilled. It is self evident that 
these formulae do not depend on the type of residual interaction, and that analogous 
expressions are obtained when the parameter g pp is used (with g pp, s for t's). 



The common behavior of the (3(3 moments for all nuclei, together with the simi- 
larity between the SMM and the full calculations for 48 Ca and 100 Mo (shown in Figs. 
|I] and |2|, respectively), suggests to go a step further and try to express the exact 
calculations within the framework of Eqs. ([)]) and (0). At a first glance this seems 
a difficult task, because: (i) the SMM does not include the effect of the spin-orbit 
splitting, which plays a very important role in the /3/3-decay through the dynamical 
breaking of the SU(4) symmetry, and (ii) the full calculations involve a rather large 
configuration space (of the order of 50 basis vectors). 

The parameters to, t\, and ti that fit the f3f3 moments displayed in Fig. [I] are 
listed in Table |2|, together with moments M.2u, ■Mo u (J 7T = 1 + ), and A4 0u (J + ^ 1 + ) 
for t = 0. The reliability of formulae @ and fllPf ) is surprising, to the extent that it is 
not possible to distinguish visually the exact curves from the fitted ones. It is worth 
noting that this situation persists even within the number projected QRPA ||. Why 
the exact calculations can be accounted for by Eqs. (0) and (|T0[)? I do not know a 
fully convinced answer. Yet, let me note that for a n dimensional configuration space, 
M.2u can always be expressed by the ratio of the polynomials of degrees 2n — 1 and 
2n in G(l+) §, i.e., 

Thus the above results seem to indicate that cancellations of the type (a) and (b) 
are likely to be operative to all orders, and that linear terms in G(l + ) are again the 
dominant ones. General expressions for M.q u , analogous to fllTD, are not known, but 
some cancellation must be taking place in these as well. 

3. Restoration of the isospin and SU(4) symmetries 

An important question in the QRPA calculations is, how to fix g pp or tl Several 
attempts have been made to calibrate g pp using the experimental data for individual 
GT positron decays. The weak point of this procedure is that the distribution of 
the f3 + strength among low-lying states in odd-odd nuclei is certainly affected by 
the charge- conserving vibrations, which are not included in the QRPA. For example, 
the single beta transitions 100 Tc — > 100 Mo and 100 Tc — > 100 Ru have been discussed 
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recently in the standard QRPA [10], where the 100 Tc states are described as pure 



pn-quasiparticle excitations, while the suggested wave function for the ground state 
in 101 Mo is (cf. ref. |TT|) is only »i 35% of quasiparticle nature: 

| 1/2+) = 0.59 | si /2 , 00) - 0.57 | s 1/2 , 20) + 0.32 | s 1/2 , 40) 
- 0.26 | 4/2, 22) - 0.26 | d 1/2 , 42) - 0.21 | g 7/2 , 24) 

(In the basis state | j, NI) the quasiparticle j and the N bosons of angular momentum 
I are coupled to the total spin 1/2.) 

We gauge t by resorting to the restoration of the Wigner SU(4) symmetry [[|. 
Unlike the method mentioned above, this method involves the total GT strength, 
which dependent of the charge-conserving vibrations only very weakly. We are aware, 
however, that the SU(4) symmetry is badly broken in medium and heavy nuclei, and 
therefore before proceeding, it is necessary to specify what we mean by reconstruction 
of this symmetry. 

For a system with N ^ Z, the isospin and spin-isospin symmetries are violated 
in the mean field approximation, even if the nuclear hamiltonian commutes with the 
corresponding excitation operators (5^ (t t and ctt t ). But, we know that when a non- 
dynamical violation occurs in the BCS-Hartree Fock (BCS-HF) solution, the QRPA 
induced GSC can be invoked to restore the symmetry. There are subtleties involved 
in the restoration mechanism: the GSC are not put in evidence explicitly, but only 
implicitly via their effects on the one-body moments (3^ between the ground state and 
the excited states. Besides, for the F excitations and when the isospin non-conserving 
forces are absent, a self-consistent inclusion of the GSC leads to the following: 

1) all the f3~ strength is concentrated in the collective state, and 

2) the 13 + spectrum, which in QRPA can be viewed as an extension of the j3~ spectrum 
to negative energies, is totally quenched. 

The self-consitency is only attained when the same S — 0, T — 1 interaction coupling 
strengths are used in the pairing and PP channels, i.e., when v pair = v pp , and the 
extent to which the above conditions are fulfilled may be taken as a measure of the 
isospin symmetry restoration. In Fig. [| is shown the behavior F strength (3 + as a 
function of the parameter s = v pp /v patr . 



Besides being spontaneously broken by the HF-BCS approximation, the SU(4) 
symmetry is also dynamically broken by the spin-orbit field and the supermultiplet 
destroying residual interactions. But, the last two effects have a tendency to cancel 
each other. In fact, within the TDA the energy differences between the GT and F 



resonances can be expressed as |12 
E GT — E F 



MeV, (12) 



2A 

where Ai s « 20A~ 1 ^ 3 is the mean spin-orbit splitting and v% h and vf h are, respectively, 
the singlet and the triplet coupling constants in the PH channel. As Vf h > vf 1 
the residual interaction displaces the GT resonance towards the IAS with increasing 
N — Z. What is more, the energetics of the GT resonances are nicely reproduced by 
(see Fig. ||) 

E GT — Ep = (26A- 1/3 - 18.5 N ~ - \ MeV, (13) 

which has the same mass and neutron excess dependence as fljjP - Briefly, the exper- 
imental data show that the SU(4) symmetry destroyed by the mean field is partially 
restored by the residual interaction. Q The GSC are likely to alter Eq. ( |T2"D very 
little. But, within the QRPA the ar + transition strength is strongly quenched and 
the GT resonance is somewhat narrowed, as compared with the TDA results. As 
such the global effect of the pn residual interaction on the GT strengths (3^ (o~t t ) 
is qualitatively similar to the corresponding effect on the F strengths /5 ± (t t ), in 
the sense that the conditions 1) and 2) are approximately fulfilled, and we say that 
the SU(4) symmetry is partially restored. It seems reasonable then to assume that 
the maximal restoration is achieved for the value of t where the GT strength (3 + is 
minimum, and this is the way how we fix the parameter t. 
4. Results for Ai eu and < m v > 

From the results displayed in Table [| it can be said that with t = t sym the 
calculated A4 ei , moment for 48 C a does not contradict the experimental limit and 
that the 2v measurement in 82 Se is well accounted for by the theory. On the other 
4 Neither in 208 Pb, where the GT strength is located at the energy of the IAS, the SU(4) symmetry 
is totally restored, as indicated by the resonance width of ~ AMeV. 
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hand, the calculated 2v matrix elements turn out to be too small for Ge and Mo 
and too large for 128 Te and 130 Te (in both the cases by a factor of ~ 3). Yet, one 
should bear in mind that: i) the calculated values of M. &u vary rather abruptly near 
t = t sym and therefore it is possible to account for the .M ej , in all the cases with a 
comparatively small variation (< 10%) of t, and ii) the minimum value of the GT (3 + 
strength critically depends on the spin-orbit splitting over which we still do not have 
a complete control. 

Besides the issue of the procedure adopted for fixing the particle-particle strength 
parameter within the QRPA, there are some additional problems in calculations of 
the matrix element M.2v, as yet not fully understood. They are related with the 
type of force, choice of the single particle spectra, treatment of the difference between 
the initial and final nuclei, etc. All these things are to some extent uncertain and 
therefore it is open to question whether it is possible, at present, to obtain a more 
reliable theoretical estimate for the 2v half lives that the one reported here. 

The upper limits for the effective neutrino mass < m v >, obtained from the 
measured 0z/ half-lives and the calculated matrix elements are shown in Table f|, 
where also are presented the results obtained by other groups. The difference in a 
factor of about 2 — 3 between both: i) the results obtained by the Pasadena group and 
the groups of Tubingen and Heidelberg for 76 Ge and 82 Se nuclei, and ii) the previous 
and present calculations for 100 Mo, 128 Te and 130 Te nuclei, is just a reflection of the 
unavoidable uncertainty of the QRPA calculations, and it is difficult to assess which 
one is " better" and which is " worse" . 
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Figure Captions 



Figure 1: Calculated matrix elements M.2v (in units of [MeF] -1 ), the Ov moments for J 71 " = 1+ 
(Mou(J 7! = 1 + )) and total moments Moi> as a function of the particle-particle S = 1, T = coupling- 
constant t. 

Figure 2: Exact (solid lines) and SMM (dashed lines) matrix elements M.2v (in units of [MeF] -1 ), 
as a function of the coupling constant t/t$. to is the value of t for which M.2i> is null. 

Figure 3: Calculated double beta decay matrix elements yV<2i/ (in units of [MeK] -1 ) for 76 Ge, 
as a function of t. Solid and dotted curves correspond to the projected (PQRPA) and unprojected 
(QRPA) results, respectively. 

Figure 4: Fermi and Gamow- Teller transition strengths f3 + for the nuclei 48 Ca, 76 Ge, 82 Se, 
100 Mo, 128 Te and 130 Te, as a function of particle-particle couplings s (S = 0, T = 1) and t (S = 1, 
T = 0) respectively. 

Figure 5: Plot of -Egt - E F versus (N — Z)/A When the experimental results overlap (for 
90,92 Zr and 208 Pb) we displace them slightly with resect to the correct value of (N — Z)/A for the 
sake of clarity. The values calculated by Eq. ( |l2| ) are indicated by full circles. 
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Tables 



Table 1: The Ml(J + ) moments and the factors G(J + )/oj° within the single mode model for 
48 Ca and 100 Mo. 



48 Ca 100 Mo 
J -Ml -G(J+)/u° -Ml -G(J + )/c° 



1.0159 



1 

2 
3 
4 
5 



1.3439 
0.1573 
0.2143 
0.0446 
0.1081 



6 0.0122 -f^.s 

7 0.0988 ±£1 



11+ 

21 L 

21 6 

19 + 
77 1 

77 * 

235 
1001 

25 
429 

175 
429 



t 



2.8316 
0.2741 
0.2086 
0.1263 
0.0585 
0.0842 
0.0177 
0.0925 



20 



t 



27 



77 \ 



130 , 

693'' 



1001 ^ ^ 81 b ' 
rt 



1001 V 1 81 
124 . 
1287 1 

143^ ^ 27 b > 
190 i. 
3861 

2431 V 1 ~ 9 a ) 



Table 2: The coefhcients to, t\, and £2 and the matrix elements M.2v, -Mov{J w = and 
A4ov(J n 7^ 1 + ) for t = 0, in the parametrization of the 2z/ and Of /3/3 moments. The matrix 
elements M.2v are given in units of [MeV] -1 . The values of the PP coupling strengths, which lead 
to maximal restoration of the SU(4) symmetry (t = t sym ), are shown in the last row. 





48 Ca 


76 Ge 


S2 Se 


100 Mo 


128 Te 


130 Te 


-M 2 u 


0.173 


0.308 


0.321 


0.451 


0.381 


0.331 


to 


1.394 


1.161 


1.206 


1.469 


1.265 


1.261 


tl 


1.754 


1.680 


1.691 


1.649 


2.131 


2.268 


-M 0v (r = i + ) 


1.506 


4.242 


4.179 


5.015 


4.599 


4.182 


-M v(J* ^ 1 + ) 


1.501 


6.924 


7.495 


9.762 


7.997 


7.486 


to 


1.227 


1.155 


1.141 


1.372 


1.377 


1.407 


tl 


1.768 


1.741 


1.764 


1.711 


2.236 


2.345 


t 2 


12.82 


13.23 


12.14 


6.527 


13.39 


11.08 


tsym 


1.50 


= 1.25 


= 1.30 


1.50 


= 1.40 


= 1.40 
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Table 3: Experimental and calculated 2v moments for t = t sym (in units of [MeV] 1 ). 





48 Ca 


76 Ge 


82 Se 


W0 Mo 


128 Te 


130 Te 


\M 2 u\ exp 
Mt 


< 0.081 
0.091 


o.28o±8:82S 

0.100 


n 1 41 +0.004 

0.121 


0.2941°;°! 
0.102 


0.0381H1 

0.118 


0.0271°;°! 
0.096 



Table 4: Upper bounds on the effective neutrino mass < m v > (in eV) obtained from the QRPA 
calculations of the nuclear matrix elements. For the sake of comparison, in all the cases the same 
experimental data, as well the same effective axial vector coupling constant (<?a = — gv) have been 
used. 
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